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Abstract. We consider the finite-size scaling of equilibrium droplet shapes for fluid 
adsorption (at bulk two-phase co-existence) on heterogeneous substrates and also in 
wedge geometries in which only a finite domain of the substrate is completely 
wet. For three-dimensional systems with short-ranged forces we use renormalization 
group ideas to establish that both the shape of the droplet height and the height- 
height correlations can be understood from the conformal invariance of an appropriate 
operator. This allows us to predict the explicit scaling form of the droplet height for 
a number of difTerent domain shapes. For systems with long-ranged forces conformal 
invariance is not obeyed but the droplet shape is still shown to exhibit strong scaling 
behaviour. We argue that droplet formation in heterogeneous wedge geometries also 
shows a number of different scaling regimes depending on the range of the forces. 
The conformal invariance of the wedge droplet shape for short-ranged forces is shown 
explicitly. 



PACS numbers: 68.45. Gd, 68.35.Rh 
1. Introduction 

There is currently mucli experimental and theoretical interest in the behaviour of fluid 
adsorption on structured substrates that are chemically heterogeneous [1-5] or that have 
a specific non-planar shape [6-10]. It is now appreciated that breaking the translational 
invariance along the substrate can lead to a number of highly interesting effects. These 
include surface phases that do not obey Youngs contact angle equation [1-3], surface 
condensation-like phase transitions [4-6], non-trivial responses of the interface shape 
to substrate corrugation 0, geometry dependent critical exponents for the growth of 
the adsorption at complete wetting [^ and even changes to the roughness exponent 
describing the fluctuations of the interface at filling transitions compared to wetting 
|l^. In the present work we are concerned with the equilibrium shape of finite-size 
droplets on planar substrates and also in non-planar wedge geometries. We argue that 
within the grand canonical ensemble (without volume constraint) the droplet shape is 
characterised by universal scaling-like functions analogous to universal order-parameter 



profiles associated with finite-size effects in critical magnets and (binary) fluids |]TT], [T^. 
In particular for systems with short-ranged forces we show that the universal droplet 
shape may be understood from the conformal invariance of an appropriately defined 
surface operator. Indeed for such systems the conformal invariance may be proved 
explicitly using a straightforward renormalization group analysis. In this way we are able 
to make predictions for the explicit droplet shape on a number of pertinent heterogeneous 
substrates including those with striped and circular domains. We beheve that our 
predictions can be tested in future numerical simulation studies of Ising models and 
may be open to experimental verification. 

Our article is arranged as follows: In §2 we discuss the case of droplet formation 
on planar substrates in (bulk) three-dimensional (c? = 3) systems (i.e. a two-dimensional 
interface) concentrating on the case of systems with short-ranged forces (for which they 
may well be experimental systems) where the conformal invariance of the droplet shape 
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and its correlations may be proved directly and used to understand the finite-size scaling 
behaviour. We finish this section with a discussion of long-ranged forces (in d = 3) 
where the droplet shape is characterised by alternative scaling functions and also a brief 
treatment of droplet formation in bulk dimension d = 2 where the result may also be 
related to the restrictions imposed by conformal invariance. In §3 we identify the scaling 
functions describing the shape of droplets in heterogeneous wedge geometries exactly at 
the filling transition temperature for systems with both short-ranged and long-ranged 
dispersion forces. Again we find that the restrictions imposed by conformal invariance 
have relevance to droplet formation for short-ranged forces. 



2. Droplets on planar substrates 

To begin we discuss the case of droplet shapes on planar heterogeneous substrates. 
Consider a planar substrate in contact with a vapour phase which is at bulk liquid- vapour 
coexistence corresponding to saturation chemical potential /x = fisat- The substrate is 
chemically heterogeneous and consists of a finite-size simply connected domain of 
substrate type A (which is assumed to extend over a region much larger than the bulk 
correlation length in all directions) which would be completely wet by the liquid phase 
if the area was infinite i.e. the local contact angle 6a = 0. Outside of the domain the 
substrate is of type B which is only partially wet by the liquid phase at this temperature 
so that 9b >0, see Fig.|l|. These conditions induce the condensation of a large finite-size 
droplet over the domain A^ with the local height of the interface pinned to a finite value 
at the edges. An extreme example of this is if substrate B has contact angle 6b = vr. 
This situation is realisable in experiments with water in which a is hydrophilic and B is 
hydrophobic say. Alternatively for Ising-like magnets one may envisage fixing the surface 
spins lying within A^ to be -|-1 and the spins outside to —1. As the area of the domain is 
increased so the equilibrium height of the droplet grows and becomes macroscopic. We 
are interested in the shape of the droplet for large domain sizes as this asymptotic limit 
is approached. In this limit the details of the wetting properties of B will not matter 
provided they serve to pin the interface to a finite height at the domain boundary. 
Similarly we do not have to make any assumptions regarding the order of any wetting 
transition occurring for the substrate. The central question that we address here is how 
the shape of the domain Aa bounding the drop affects the equilibrium droplet shape. 
Note that there is no volume constraint imposed on the amount of liquid adsorbed on the 
domain so the equilibrium shape of the interface necessarily has the same translational 
invariance as the underlying substrate binding potential (see 0] for effects induced 
by volume constraint or equivalently by over-saturating the bulk vapour phase). The 
simplest possible starting point for discussing this problem which includes the effect of 
interfacial fluctuations and microscopic interactions is the interface Hamiltonian model 

H[e] = /dx{|(V£)^ + l^(£;x)} (1) 
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where £(x) denotes the local interface height at position x= (x, y) along the substrate, S 
is the interfacial stiffness (tension) of the unbinding fluid interface and W{£; x) denotes 
a local position dependent binding potential. As the position x changes from outside 
to inside the domain so the local potential changes from being of type WB^i) to VFa(^) 
as appropriate to the respective homogeneous B and A substrates. The cross-over in 
W{i; x) can be anticipated to occur over a finite distance A^^ @) which itself depends 
on the range of the forces in the system. For short-ranged forces we anticipate that b 
is of order the bulk correlation length. This model also assumes that the interfacial 
gradient is small otherwise the first term in (|l]) has to be replaced with the usual 
drumhead expression for the interfacial area. We will only be interested in the scaling 
properties of the droplet shape however which emerge when the gradient is small. For 
the model (|l|) it is straightforward to derive an exact Ward identity for the equilibrium 
interface position. Following the analysis of Mikheev and Weeks [|13| (see also fl^) it is 
straightforward to show that 

= {dW{£;x)/d£) (2) 

which relates the local effective force to the local interface curvature. Here (■) denotes 
the usual ensemble average. Thus provided we restrict our attention to the scaling limit 
of domain area — > oo and positions x which are far from the domain edge then the 
shape of the droplet is described by the solution of 

^V'{i) = {W'Ai)y, xgAa (3) 

which serves as the starting point for our analysis. This scaling limit is analogous to that 
encountered in studies of universal order-parameter profiles in confined critical systems 
in, O. To proceed we use the very well developed renormalization group (RG) theory 



of wetting [15-19] and replace the RHS of (H) with the derivative of an effective potential 
which is constructed on the basis of the linearised RG. In other words we rewrite (|^) as 

SV^(£) = iy,V/(£)); xgAa (4) 

where the effective potential is derived by coarse-graining the bare potential over the 
local roughness of the interfacial fluctuations. The net result (in dimension d = 3) is 
[15-19] 

/■oo -{^-t)2/2ei 
Weffie)= dtw{t) — - (5) 



where ^± = y is the local effective roughness and is a transverse correlation 

length determined self-consistently by C,^'^ oc W^'jrj-{i). This then is the desired field 
equation for the mean droplet height, the precise form of which depends on the range 
of the forces. Note that the existence of the effective potential implies that there exists 
a local free-energy functional F[i] for the mean interfacial height (and its correlations) 
which is of the same form as (|T|) (at least within the desired region x G A^) but with 
WA^i) replaced with Weff{i)- The existence of a local free-energy functional for the 
droplet problem is directly related to the fact that the exponent analogous to the bulk 
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correlation function exponent rj is zero for wetting transitions. This has important 
simplifying consequences for the structure of correlations in the droplet which we will 
discuss later. 

To continue we consider the cases of short-ranged and long-ranged forces separately. 



2.1. Droplet shapes for short-ranged forces 

For systems with short-ranged forces such as Ising and Landau-Ginzburg- Wilson models 
and also some experimentally accessible systems including polymers ||20[, some fluid 



mixtures |^ and superconductors p2[ it is well understood [15-19] that the bare 
(unrenormahsed) potential decays as WOi(^) ~ Ase~'^^ with Ag an effective (positive) 
Hamaker constant and k = 1/^b is the inverse bulk correlation length We 
emphasise here that because the droplet shape reflects finite-size effects occurring at 
complete wetting we only need account for the leading order term in the decay of the 
binding potential. The renormalization of this is very well understood and gives rise to 
an effective potential which depends on the value of the wetting parameter cu = . 
For values of < 2 as is pertinent to fluids and Ising systems there is no need 

to consider the "hard wall" restriction / >0 (which can only be treated approximately 
by the linear RG method) and the effective potential is Ase~'^^ where k = n/{l + uj/2) 
is a renormalised inverse lengthscale. The field equation for three-dimensional droplets 
with short-ranged forces is therefore 

SV'(£) = -Aske-^'^'^ (6) 

which we wish to solve for different shapes of the domain A^. The central observation 
of this paper is that the above equation is conformally invariant with e"'^^^^ playing 
the role of a primary operator with scaling dimension Xb = 2 p5|. The proof of this 



is straightforward, but tedious and is presented in the appendix. Specifically under a 
conformal map z^w, where w = u + iv is an arbitrary analytic function of the complex 
variable z = x + iy, the field equation (^ is invariant provided we identify 

To exploit this result consider the simplest example of a structured substrate consisting 
of two semi- infinite planes of type A, B which meet at the line |/ = (See Fig.^. From 
(^) it follows that e"''^^^ shows the scaling behaviour 

so that fi;(£)oo/2~2 ln(?/C~^/^). Here C = 2Tj/ Ak^ is a non-universal constant. Following 
the analysis of Burkhardt and Eisenreigler we now use some well known analytic 



functions to conformally map the semi-infinite system onto a number of pertinent finite- 
size domains and use (^) to identify the appropriate droplet shape. For example the 
analytic functions w = -Inz and w = - cosh~^ z map the half plane onto the infinitely 
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long strip {— oo<n<oo} U {0<v<L} and half-strip {■u>0} U {0<f <L} respectively. 
In this way we obtain the droplet shapes as 

^C'f-sin^V' (9) 



,7r L 

for < y < L and 



'^{^Z half strip _ ^ 



L Tcx\ f L , Try 



sinh — + — sin 

vr L / \TT L 



(10) 



for < y < L and x > 0. Note that in writing these expressions we have changed our 
co-ordinate notation back to {x,y) rather than persisting with {u,v) (see Fig.|^). 

For circular domains of radius R the special conformal group can be used to obtain 
the desired droplet shape 



-I^Wdisc — (J 



I 



r \ 



n ^ -2 



(11) 



MJ 

where r denotes the distance from the centre of the disc. Note that the special conformal 
group can be used to map hyperplanes onto hyperspheres in arbitrary dimension so 
that we may anticipate that a result analogous to (p!TD may exist in other dimensions 
in particular if the substrate (or unbinding interface) is effectively one dimensional. We 
shall return to this point later. 

Equations (^ and (0) are our main predictions for droplet shapes in systems with 
short-ranged forces. We believe these to be the most relevant for future numerical 
and simulation studies and we finish this section with a discussion of them. Our first 
remark is that the droplet-shape shows the same kind of finite-size scaling behaviour as 



predicted for order-parameter profiles in confined critical systems |0, 12 . Consider for 



example the result (|^) for the strip domain. For large (but finite) distances from the 
droplet edge at the y = line, the local height shows the " distant wall correction" 

2 2 

mstrip = Hi)oo/2-Y{l) +••■ (12) 

for y/L'^1. The power of L here reflects the two-dimensional character of the interface 
|11[|. Secondly our results (||) and (|Tl|) for the droplet shape can be conveniently recast 
in scaling form by re-writing them in terms of the maximum height i^"^^ of the droplet 
occurring at the middle of the strip, disc domains. Thus for the strip we predict 

-4::U = (2 + ^)lnsin^ (13) 



/strip '•-strip J 

whilst for the disc 



^ {{^Usc - iSc) = (2 + 0.) In f 1 - ^) . (14) 



Again we emphasise that the scaling function ( |T3| ) is valid for y oo, L ^ oo with y/L 
arbitrary (and similarly for (p!^). It is these expressions that we believe would be easiest 
to test in future simulation or experimental studies. Before we discuss the possibility of 
this we mention here that in contrast to studies of order-parameter profiles in confined 
critical magnets and fluids where the absolute value of the singular contribution to the 
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mid-point magnetisation (say) is very difficult to measure (because it vanishes as L — oo 
or i?-^oo ) the scahng of the droplet height is an obvious experimental observable. In 
particular the above results make a precise prediction for the difference in the maximum 
(mid-point) droplet height for the strip and disc domains. Specifically the difference in 
height of a droplet on an infinite strip of width L and that on a disc of identical diameter 
satisfies 

'^(4::l-4r2) = (2 + ^)ln^ (15) 

which depends only on the value of the wetting parameter. It is remarkable how small 
this height difference is. In mean-field approximation which ignores fluctuation effects 
and has an effective value of 07 = the numerical value of the RHS of (piSf) is only 0.4831, 
i.e. the droplet heights are different by less than half a bulk correlation length. Close to 
the bulk critical point where u approaches a universal value tUc — 0.77 this number 
increases to 0.67 but again the difference between the two heights is nearly negligible. 

One of the ffist tests of these predictions would be to numerically study a mean-field 
Landau-like model of adsorption on heterogeneous walls. The influence of fluctuations 
on the droplet shape is rather small and the above results are valid even in the mean- 
field limit of to' = 0. Beyond mean- field level it would also be straightforward to consider 
Ising model simulation studies similar to those performed for wetting at homogeneous 
walls We note that for this case it may well be that the appropriate value of 

k is altered by the coupling of interfacial fluctuations p9|, pOl. This is certainly needed 



to quantitatively explain simulation studies of (other) finite-size effects at complete 
wetting One obvious experimental constraint on these predictions is the role of 

gravity which both thins the height and dampens the the fluctuations of the liquid- 
vapour interface. However this effect will only become important if the strip width 
(say) is larger than the capillary length for gravity damped interfacial modes. Using the 
value of the surface tension of water at room temperature this implies that gravity may 
be neglected for strip widths or disc radii of less than a few micrometers. 



2.2. Height fluctuations and correlations in the droplet 

Next we turn our attention to the scaling properties of the height-height correlation 
function S'(xi,X2) = (£(xi)£(x2)) — (£(xi)) (£(x2)). There are a number of ways of 
calculating this quantity but following the treatment above we use the effective potential 
(^) to derive a renormalised Ornstein-Zernike (OZ) equation which for short-range forces 
reads 

(-SV,^ + Ak^e-'^') 5(xi, X2) = 5(X2 - xi). (16) 

Here we point out the invariant properties of this equation as well its explicit solution for 
the semi-infinite and strip geometries. On making use of the covariant relation (^ it is 
straightforward to show that this equation is invariant under a conformal map z^w{z) 
with the correlation function transforming simply as 

S'(ui,U2) = 5'(xi,X2) (17) 
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where we have adopted the notation ui = {ui,Vi). The absence of an exphcit scahng 
operator term in this resuh reflects the marginal logarithmic growth of the droplet 
height away from the boundary for systems with short-ranged heterogeneous forces. 
The transformation law for the height correlations is equivalent to the statement that 
the (connected) two-point correlator of the primary operator e~'^^ transforms as 

^g-K^(u,)g-K^(u.)^^, = |w'(;2i)|-V(-S2)r'(e-"'("^)e-"'^("^))G (18) 

where the subscripts G' and G refer to ensemble averages with respect to the different 
geometries under the conformal transformation. That ( p!7D and (|1^) are equivalent again 
owes itself to the existence of a renormalised local binding potential which accounts 
for fluctuation effects. Having established that the height-height correlation function 
transforms conformally it immediately follows from the analysis of Cardy that on 



the semi-infinite planar surface (Fig.^ the two-point function depends only on a single 
variable such that 

The scaling function is then determined by substituting this result into the OZ equation 
(|T6|) with e~'^^ given by (0). We find that <l>(t) satisfies 

(t^ - l)$" + 2t$'-2<l> = (20) 

from which we can easily determine the important asymptotic behaviour. At short 
distances xi X2 corresponding to t ^ 1 it follows that $ ~ ln(t — 1) implying that 
the local roughness ^± ~ ^Iny precisely as expected for interfacial fluctuations at the 
marginal dimension d = 3. On the other hand in the limit of large transverse separations 
|x2 — xi| — >oo corresponding to t— >-cx) it follows form ([l3|) that which identifies 



the surface scaling exponent t]\\ =4 |3^, It is this limit of the scaling function which 
is required in order to understand the asymptotic decay of correlations in the infinite 
finite width strip geometry. Using the mapping discussed above it immediately follows 
from ([T2|) that the height-height correlation function decays as 

5.,„,(xi,X2) ^ sin^ (^) sin^ (^) e"^-!---!/^ (21) 

for — Xi\^ L. In writing this expression we have again chosen to use the original 
co-ordinate notation rather than {u,v). The above expression identifies the universal 
finite-size droplet correlation length = L/2tt. This identification is consistent with the 
amplitude-exponent relation L/^y =2/7r?7|| |Tl[] . 

We note here that the same result for 5'(xi,X2) in the strip follows from explicit 
solution of the OZ equation (|TTp using the derived expression for the droplet height (^. 
This can be written clS cl spectral expansion 



5.t„p(xi,X2) = E C(2/i)^n(l/2)e-^"l^^-^^l (22) 

ra=0 

with eigenvalues and eigenvectors determined from solution of 

-< + 2(^sin^) %n = i?^^n (23) 
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subject to ipn{0)=4'niL) = and the orthonormality condition 



L 



2S / dyrniy)My) = SnmE-\ (24) 







The solution to (p3D is found using standard methods. For the eigenvalue spectrum we 
find En = {n + 2)'k / L (with n = 0, 1, 2, . . .) whilst the eigenvectors are 



n n 1 2 TTX 
— , 2 H — , -; cos 
2 2 2 



(25) 



1 n 5 n 3 2"^^ 

2 ~ 2 ' 2 + 2 ' 2' T 



(26) 



for = 0, 2, 4, . . . and 

T/'nlx) = iV„ sm ( — I cos —F 

for n = 1, 3, 5, . . .. Here F[a,/3,7;x] denotes the usual hypergeometric function. We 
note here that the fact that the exact result for S'(xi, xi) in the strip geometry obtained 
by conformally mapping the the semi-infinite correlation function also follows from the 
solution to a local field equation (|16D is intimately related to the existence of a local 
effective free-energy functional for complete wetting which in turn reflects the fact that 
there is no exponent analogous to the rj occurring for bulk critical systems. In studies of 
flnite-size effects in strip (and other) geometries occurring at bulk criticality (for which 
77 7^ 0) local fleld equations for the pair correlation function derived from approximate 
local functional models do not capture all the necessary scaling features required 
by conformal invariance (although they are surprisingly accurate and highly reliable for 
certain quantities ||35[| ). 

As a flnal remark concerning height-height correlations in the strip geometry we 
note that the existence of a renormalised local binding potential also implies that the 
correlations in the strip satisfy the algebra of correlations discussed by Parry and Swain 
36| in the somewhat different context of magnetisation functionals. Deflning the Fourier 



transform of the correlation function S{yi,y2', q) by 

^(z/i,l/2;q) = / dx2ie-'^(^-''^)S(xi,X2) (27) 

and in particular denoting the zeroth moment So{yi, 1/2) = S{yi, 2/2; 0), it follows from the 
local character of the renormalised effective free-energy functional that for any ordered 
points 0< yi< y2<y3< L, the correlations satisfy 

S{yi,y2;q)S{y2,y3;q) = 5(?/i,2/3;q)^(2/2,2/2;q) (28) 



and 



So{yi,y3) 5*0(1/1, 1/2) \ / Soiyi^ys) So{y2,y3) 



.('{yiWiys) nyiy'{y2)) V{yiy'{y^) ^'(y2)^'(i/3) 

Soiyi^yi) 5*0(1/1, 1/2) \ / So{y3,y3) 5*0(1/2, ys) 



(29) 



V'iyi)i'{yi) i'{yiy'{y2)J V'{y3y'{y3) i'{y2)i'{y3) J 

where i'{y) denotes the derivative of the equilibrium height proflle given by (^. These 
relations severely restrict the possible form of 5o(?/i,y2)- In particular they imply that 
the dimensionless quantity 
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where AF{L) denotes the excess or finite-size contribution to the droplet free-energy, 
satisfies 

(^ivi, y2) = (i - \/i - ± \/i - fT(y2,z/2)) ; 1/1 < § (31) 

and the ± signs refer to whether yi and y2 are on the same/opposite side of the mid-hne 
at y = L/2. We emphasise here that the results (p8D-(pl|) are valid for any choice of 
effective binding potential and are not specific to the exponential form for systems 
with short-ranged forces. The central prediction (|3T| ) of the correlation function algebra 
relates the correlation function structure within the droplet to the excess free-energy 
AF{L). For the droplet with short-ranged forces AF{L) is readily calculated using the 
appropriate effective free-energy functional. We find 

AF(L) = —{2 + luYL'\ (32) 

Again the width dependence of this quantity is precisely in accord with dimensional 
(hyperscaling) ideas for confinement in critical systems [|12|. Note that the (Casimir) 



amplitude here is strongly non-universal and depends on the value of u. This is in 
contrast to finite-size effects at bulk criticality (in Ising systems for example) where the 
amplitude of AF{L) is universal and related to the central charge characterising 



the bulk universality class. The derivative of the free-energy determines the force of 
solvation between the edges of the strip which might be measurable in force-balance 
experiments if the width L of the droplet can be controlled externally. 

2.3. Droplets shapes and correlations with long-ranged forces 

In the last section we explicitly showed and exploited the conformal invariance of the 
droplet shape and its correlations for systems with short-ranged forces. Whilst some 
examples of such systems may well be found in the laboratory most solid-fluid and 
fluid-fluid systems will de dominated by dispersion forces for which the bare binding 
potential (of the completely wet A region) decays like Wa ~ where again Ai is an 
effective (positive) Hamaker constant. It is well understood [37| that this potential is 



not changed by fluctuations in c? = 3 so that the fleld equation for the droplet shape is 
simply 

EV^i) = -2Ai{i)-\ (33) 

This equation is certainly not invariant under the same kind of conformal mapping 
as valid for systems with short-ranged forces consistent with general expectations 
of conditions for conformal invariance. In order to investigate the possible scaling 
properties of different flnite-size droplet shapes it is therefore necessary to solve (^) 
explicitly for different domains. Here we concentrate on the inflnite strip and disc 
domains only. For the strip geometry the differential equation is trivially solved and 
shows that the equilibrium height proflle has a simple scaling form 

{i)str^, = 24"Wf fl-y); 0<y<L (34) 
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where i^^lp oc L^/^ is the maximum (mid-) height of the drop. Thus whilst there is no 
conformal invariance for the droplet with dispersion forces there is still a strong degree 
of scaling in this system. Indeed the scaling of the height profile for long ranged forces 
is stronger than that predicted for short ranged forces as can be seen from comparison 
with ([T3|). Notice that near the edges of the droplet the interface height behaves as 



^strip^y^^"^ (for y/L^O say) where the power law here can be identified with the ratio 
of complete wetting exponents Pco/ ^\\ where Pco is the adsorption critical exponent and 
z/|| is the transverse correlation length critical exponents. For systems with dispersion 
forces recall that /5co = 1/3 and i^y =2/3 consistent with the explicit shape dependence 
give above (see for example [|T^ ). 

Correlations along the droplet show a similar scaling structure. From the OZ 
equation 

' -SV,2 + '5(Xi,X2) = <5(X2 - Xi) (35) 

it follows from the spectral analysis given earlier that 5'(xi,X2) decays asymptotically 

as 

5(xi,X2) ~ ^o(2/i)^o(y2)e-^»l^'^-^^l; ^ oo, (36) 

where tpo and Eq are found from numerical solution of the eigenvalue problem 

3L^ 

4^2(^2,^2) ^o = i^g^o (37) 

with ipo{^) = "^oiL) = 0. In this way we find that the finite-size correlation length 
characterising height fluctuations along the droplet satisfies L/^\\ ~ 5.25 which is 
universal in the sense that it does not depend on the value of the surface tension or 
Hamaker constant. Interestingly the ratio is smaller than that derived earlier for short- 
ranged forces in keeping with general expectations about the influence of long-ranged 
forces on the range of correlations. 

For droplet shapes on circular domains of radius R it follows from (|33D that the 
equilibrium height profile also has a scaling form {i)disc = ^disc-^i^/-^) where £^^^^oci?^/^ 
is the mid-point height and A(^) is a scaling function satisfying A(0) = 1 and A(l) = 
which can be determined numerically. Following our earlier treatment of droplets with 
short-ranged forces it is interesting to compare the maximum droplet height for a circular 
domain of diameter 2R = L with the corresponding droplet height for a strip geometry. 
Because there is a stronger manifestation of scaling for systems with long-ranged forces 
the comparison takes the form of a ratio £^™)p/£^^),^ 1.03 rather than a difference as for 
short-ranged forces (|15|). The influence of the domain shape on the droplet adsorption 
is therefore stronger for systems with dispersion forces even though the scaled heights 
are again very close to each-other. A comparison of the corresponding droplet shapes is 
shown in Fig.H for strips and discs with long-ranged forces. 
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2.4- Droplet shapes m 1 + 1 dimensions 

To end this section and motivate out treatment of droplet shapes in three-dimensional 
wedges we remark here that the scaling form of the droplet height for adsorption on 
a one- dimensional heterogeneous substrate also shows scaling behaviour that can be 
understand using conformal invariance. To be concrete consider a two dimensional 
semi-infinite Ising model at temperature T (less than the bulk critical value T^) and 
in zero magnetic field H = 0~ (so that the bulk phase has negative magnetisation). 
The surface spins on the y = line are all fixed to be down except along the line 
between x = and +L where they are all up. This is the two-dimensional analogue 
of the striped heterogeneous surface and the circular domain considered above. These 
boundary conditions induce the formation of a droplet of local equilibrium height i{x) 
which is pinned to the wall at the end points x = and x = L. In the limit of L^oo the 
droplet becomes macroscopic in size and we can anticipate that the shape shows scaling 
behaviour. The interfacial model for this problem is 




where the binding potential W{i) may be simply approximated by an infinite hard- 
wall potential W{i) = for i> and W{i) = oo otherwise, which simply restricts the 
interface to the / > region. The reason for this is that in d = 2 the entropic repulsion 
of the interface from the hard-wall dominates over any short-ranged interaction (for 
complete wetting problems) [17-19]. This problem has been considered by Burkhardt 
using continuum transfer-matrix techniques with the scaling behaviour of £{x) 
explicitly calculated to be 

W^)) = 2£ryf (l-f) (39) 

where i'^^ oc L^/^ denotes the mid-point droplet height. Note that this expression is 
identical to that derived above for 2 + 1 dimensional systems with dispersion forces. 
Here we make two relevant remarks concerning this transfer matrix result : 

(i) the scaling of {i{x)) in two dimensions is precisely what one would expect due to 
the entropic repulsion from the hard-wall. The droplet shape reflects interfacial finite- 
size effects in the weak-fluctuation regime (WFL) |1^, [1^ which can be quantitatively 
understood using an effective potential which decays as Wfl'^^fl^~'^- Using this form 
for the effective potential the field equation (^) reduces to 

= -2A^i(£)-3 (40) 

which is trivially integrated to recover the exact transfer matrix result (pQ]). 

(ii) the scaling of {i{x)) is also in accord with the predictions of conformal invariance 
for (non-trivial) one- dimensional field theories. In general phase transitions do not 
occur in (classical) one-dimensional systems unless either extremely long-ranged forces 
are present or the order-parameter is allowed to diverge. The study of droplet shapes 
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in 1 + 1 dimensional systems is therefore an interesting test of the ideas of conformal 
invariance in this dimension. To see this first consider the limit of L ^ oo. Then far 
from the pinned edge of the droplet the height i{x) scales as 

{i{x)) ~ Dx'^; Kx > 1 (41) 

where C, = 1/2 is the usual interfacial wandering exponent for the one- dimensional 
interface [|T^, |T^. This behaviour is consistent with the homogeneous rescaling x—>-x/b 
under- which {i{x)) transforms as 

{i{x/b)) = b-^{i{x)). (42) 

It is natural to generalise this statement to non-uniform rescalings corresponding to the 
one-dimensional limit of special conformal mappings such that 

/dr'\'^ 

Wg' = (^^j m)G, (43) 

where as before G' and G refer to ensemble averages with respect to the old and 
new geometries respectively. The allowed class of transformations combines uniform 
rescalings, translations and inversions is therefore of the form 

x' = — (44) 

c + dx ^ ^ 

with ad — bc^O (since otherwise the scale factor vanishes). This transformation maps 

the semi- infinite line to a line of finite length. Choosing a, b, c and d such that < a;' < L 

we find 

( 



{i{x')) = £('")4^ 




(45) 



which is of course identical to the transfer matrix result (^) provided we identify the 
correct value of the wandering exponent. Note that if we write this expression in terms 
of the distance x" = x' — L/2 from the centre of the strip then the position dependence 
of the scaling function is basically identical to that observed for circular domains in 
d = 3 but with a different value of ( (see (|l^) ). This is because the scaling function 
for the height of radially symmetric drops on heterogeneous substrates (with short- 
ranged forces) in arbitrary dimension can be obtained by the application of the special 
conformal mapping analogous to (||). 

To finish this section we emphasise here that the only ingredient in deriving the 
prediction ( P5| ) is the assumption of co-variance under the special conformal mapping 
(^). The finite-size droplet in 1 + 1 dimension is a specific example of this one 
dimensional prediction of conformal invariance for a problem in which ( = 1/2. In the 
next section we shall show that ( |15| ) has application for a second example of interfacial 
phenomena. 
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3. Droplet shapes in three-dimensional wedge geometries 

The calculations presented in section §2 assume that the substrate is heterogeneous 
but still planar. However recently it has been shown that chemically homogeneous but 
non-planar substrates in the shape of wedge geometries exhibit filling (or wedge- wetting 
transitions) which are characterised by very strong interfacial fiuctuations belonging to 
quite distinct universality classes and fiuctuation regimes to those occurring for wetting 
at a planar wall |T^. The purpose of this section is to investigate the scaling of droplet 
shapes in wedge geometries. Before we study these finite-size effects we briefly review 
the phenomenology of flUing transitions and the recently developed fluctuation theory 
for them based on analysis of a novel interfacial Hamiltonian . This model will form 
the starting point for our study of droplet shapes in wedges. 

3.1. Filling transitions and effective Hamiltonian theory 

Consider a wedge formed by the union of two walls at angles ±a to the horizontal. The 
wedge is supposed to be in contact with a bulk vapour phase at saturation chemical 
potential fi = fisat{T)~ and is at a temperature T less than the wetting transition 
temperature of the planar wall-vapour interface. Thermodynamic arguments [39- 
41] dictate that the wedge is completely filled with liquid at a temperature Tp satisfying 

eiTp) = a (46) 

and 9(T) denotes the contact angle of the liquid drop on the planar substrate. This 
thermodynamic prediction is confirmed by specific model calculations which include 
the infiuence of intermolecular interactions and fiuctuation effects |10|- The filling 
transition occurring as T — Tp may be first or second-order |^ corresponding to 



the discontinuous or continuous divergence of the interfacial height as measured 
from the bottom of the wedge. Importantly even substrates that exhibit first-order 
wetting transitions can exhibit continuous filling transitions [|10] so that the experimental 
observation of critical (continuous) filling is a realistic possibility. In addition to the 
divergence of £o as T — > Tp , the filling transition is characterised by the divergence 
of distinct correlation lengths along C,y ~ (Tp — T)~'^y and across ~ (Tp — T)~'^"= the 
wedge, as well as the roughness of the unbinding interface (see Fig.|) |jlO|. The mean- 
field critical behaviour has been examined in detail [|lOl from analysis of the non-planar 
version of the standard interfacial Hamiltonian (valid for small a) given by 

H[l] = I dxl^^iVif + W{i-a\x\)^ (47) 

where £(x) is the local height of the interface as measured from the plane and W{1) is 
the usual binding potential. The mean-field critical exponents for filling are determined 
only by the leading order decay term in W{i). Writing this as W{i) = —Aii~P (so that 
p = 2 corresponds to non- retarded dispersion forces) the critical exponents are given by 
1 1 11 1 

p p 2 p 4 
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Notice the remarkable universal value of the roughness critical exponent. Analysis of the 
height-height correlations at filling show that fluctuations across the wedge are highly 
localised to the filled region with the local height £(x) at position y along the wedge 
fiuctuating coherently, i.e. the effective stiffness across the filled portion of the wedge is 
infinite. In this sense fluctuations across the wedge are trivial which can be seen in the 
explicit values of the critical exponents since (3q = for all values of p. On the other 
hand correlations along the wedge are characterised by a soft mode and completely 
dominate the fluctuation behaviour. On the basis of this Parry, Rascon and Wood 



have argued that the fluctuation theory of three-dimensional fllling transitions can be 
understood using the effective one-dimensional model 

H,yi,^^jJ^[f;^V(Q) (49) 



where £o(y) =^(0, y) is the local height of the interface at the centre of the wedge. This 
model is based on the assumption that only fluctuations along the wedge contribute 
to the asymptotic critical behaviour (see Fig.^. This effective model is only presumed 
correct for the soft-mode fluctuations in the y direction which determine the asymptotic 
critical behaviour. The model is not valid for wavevectors Q > Qmax — ^Hx- At this 
scale the full 2d model ( ^7]) has to be used. Notice that the coefficient of the bending 
term is proportional to the local height of the filled region (quite unlike the stiffness 
term for planar fiuid interfaces) refiecting the fact that a local increase in the height 
also increases the width of the filled region in the x direction. The direct filling potential 
appearing in (|i9|) is found to be 



V{^^) = -{9' - a') + Apf-P + ... (50) 
a 

the minimum of which identifies the mean- field £o. This filling Hamiltonian can be 
studied exactly using transfer matrix techniques yielding a complete classification of the 
possible critical behaviour. For p < A mean-field theory is valid, whilst for p > 4 the 
critical behaviour is universal with critical exponents given by 
113 1 
/3o = ^ ; = 4 ; '^y = l " 4' 

Notice that in the fiuctuation dominated regime /3o = z/_|_ so there are regular excursions 
of the interface to the wedge bottom. The existence of two fiuctuation regimes for critical 
filling can be understood using arguments analogous to the well developed scaling theory 
of wetting [17-20]. To see this we suppose that the infiuence of fiuctuations at critical 
filling can be understood using an effective potential 

Veffih) = V{io) + VpLiio) (52) 

which accounts for fiuctuation effects by adding to the direct potential a term arising 
from the entropic repulsion form the wedge bottom. The fluctuation term can be 
estimated from the form of the bending energy contribution to (^91) as Vfl(£o) — ^o^y ^ 
and hence 

VpLiio) ^ ^o~'^^ (53) 
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where ( is the wandering exponent for filhng transitions which relates the roughness 
^± ~ and soft-mode correlation length ^y. The value ( = 1/3 can be either read 
directly from the transfer matrix results (^Tj) or equivalently from the scale invariance 
of the free filling Hamiltonian Hq = J dy^ (^) ^^'i^'^ the scale transformation y^y/b 
and i^i/b'^. The net result of this is an effective filling potential of the form 

V,ffiio) = Vi£o) + Di^' (54) 

which is minimised to find the equilibrium film thickness £o. Thus we find a mean- 
field and fluctuation dominated regime depending on whether the entropic repulsion 
or direct algebraic term Api^~^ is the dominant correction to the linear term. Unlike 
critical wetting for which there are three fluctuation regimes critical filling only has two 
since the linear term in the direct binding potential is always relevant. In this sense 
fluctuation effects at three dimensional critical filling is rather similar to those occurring 
at two-dimensional complete wetting [17-20,37]. 

3.2. Finite-size droplets in heterogeneous wedge geometries 

We are now in a position to consider the finite-size scaling of droplet shapes in 
heterogeneous wedge geometries formed by sandwiching a finite length wedge of 
substrate type A between two wedges of type B (see Fig.^). Following our earlier 
treatment of heterogeneous planar substrates we suppose that exactly at saturation 
chemical potential /i = fisat{T)~ the filling height on the B wedge is finite but that 
an infinite domain of A wedge would be completely filled. This geometry induces the 
condensation of a large droplet of liquid over the A region with the interface height 
pinned to finite values at the end points (taken to be |/ = and y = L). As before there is 
no volume constraint on the amount of liquid adsorbed so that as the maximum height 
of the droplet oo as L ^ oo and we can anticipate scaling behaviour for the 

position dependence of the equihbrium profile {£o{y)). The details of the partial filling 
of the B region are not important for the scaling structure of the droplet shape and 
there are a number of specific ways of inducing the pinning. For example if the contact 
angle of the liquid drop on the planar B substrate 6*^ > 7r/2 then the associated wedge is 
never filled (by liquid) for any temperature. If 9b <vr/2 then one has to chose substrate 
A such that 9a <9b so that the A wedge fills before the B wedge. We now distinguish 
between two distinct manifestations of scaling for the finite-size droplet shape: 

(A) Droplet shapes exactly at the filling transition temperature. 

For this case the temperature T is tuned exactly to the filling transition of the 
fluid on the (infinite) A wedge, i.e. the temperature satisfies 9a(Tp) = a. This is the 
marginal value of the temperature at which the maximum droplet height i]^""^ becomes 
macroscopic as L—>-oo. 

(B) Droplet shapes above the filling transition temperature. 

If T>Tp then the droplet height also grow to infinity as L^oo. However because 
of the relevance of the linear term in the filling binding potential (|50D we cannot assume 
that the scaling of the droplet shape is the same as that obtained exactly at T = Tp. 
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Indeed droplet shapes for T > Tp may be viewed as finite-size effects occurring at 
the complete-fiUing transition. However fluctuations here do not have the same quasi- 
one-dimensional characteristics of critical filling. Analysis of the critical exponents 
at complete filling show that and diverge with the same critical exponent as 
fi^fi~at{T) 1^, [ly]. Consequently the 1-d model (^Qf ) is not justified for this transition 
and analysis of the droplet shape in the present heterogeneous wedge problem must be 
based on the full model (|47D . This is a much more difficult problem and we will not 
consider this here. 

Our treatment of droplet shapes in the wedge parallels that of droplet shapes on 
planar substrates. The effective Hamiltonian on the heterogeneous wedge is taken to be 
the generalisation of the filling model (^) 

where V{£; y) is a position dependent direct filling potential that changes from being of 
type Va(^) to Vb(^) as one crosses the heterogeneous region. The Ward identity for this 
model is again easily derived. We find 

^ {io{y) Qy)) = (^^{i; y)^^ (W) (56) 
so that within the droplet we can write 

— (Uy) ay)) - - i^y?) = (vm- (57) 

a a 

To evaluate the RHS we again resort to the concept of an effective potential which we 
showed to yield exact results for the wetting droplet in 1 + 1 dimensions. Thus we seek 
to solve the field equation (writing (£9) = ^o, for convenience) 
9V y 

— ^o^[; + -^[,^ = K'//(^o) (58) 
a a 

using the effective potential K// evaluated at T = Tp for short-ranged and long-ranged 
forces. For long-ranged forces we will assume that p = 2 as is appropriate for realistic 
solid-fluid interfaces. This differential equation is easily integrated using standard 
techniques and we only quote the final results. 

For droplet scaling exactly at Tp and short-ranged forces we find 



(£o(?/)) = C"'4i/3 



y (^_y 



nl/3 



(59) 



L V L, 

where P^""^ (xL^^^ denotes the mid-drop height. Note that this scaling form is precisely 
in accord with the prediction of conformal invariance for one dimensional interfaces with 
wandering exponent C = 1/3. For long-ranged forces on the other hand integration of (^) 
yields a quite different form of the scaling function. We find that (^0(2/))/^™"^ = ^o(y) 
satisfies the cubic equation 

3£l + £^ = i6^fl-^) (60) 



with Thus the scaling of droplet shapes at wedge wetting is only consistent 

with conformal invariance for systems with short-ranged forces. 
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4. Conclusion 

In this paper we have investigated the scahng behaviour of finite-size droplet shapes on 
heterogeneous planar substrates and wedges. Our predictions are restricted to systems 
at two- phase bulk liquid- vapour coexistence (chemical potential n = iXsat) in the grand 
canonical ensemble, i.e. without volume constraint. The main conclusions of article are 
the following: 

(A) For three-dimensional systems with short-ranged forces the scaling of the 
droplet shape and its height-height correlations can be understood using conformal 
invariance. In fact using standard wetting renormalization group ideas the conformal 
invariance may be explicitly proved for such systems. Consequently the study of 
droplet shapes in laboratory systems with effective short-ranged forces such as polymers, 
superconductors, liquid metals and even some fluid-fluid interfaces may be viewed as a 
possible experimental testing ground for many of the ideas developed in conformal-fleld 
theory as well of course as the fluctuation theory of fluid interfaces. 

(B) For three-dimensional systems with long-ranged (van der Waals) dispersion 
forces the droplet shape still exhibits a strong degree of scaling invariance related 
behaviour even though there is no longer any conformal invariance. 

(C) For model systems of interfaces in two-dimensional bulk geometries with 
short-ranged forces the scaling shape of the droplet on the heterogeneous wall can 
be understood quantitatively using an effective potential accounting for the entropic 
repulsion of the interface from the wall and is completely in accord with the predictions 
of conformal invariance for one dimensional fleld theories. 

(D) Finite-size droplets in wedge geometries show a number of interesting scaling 
behaviours depending on whether the interfacial pinning is done exactly at or above 
the filling transition temperature. Exactly at Tp and for systems with short-ranged 
forces the prediction of conformal invariance for the shape function is again valid. For 
long-ranged forces a different scaling behaviour is observed. 

The predictions of this paper arc based entirely on the analysis of effective interfacial 
models. It would be interesting to test these ideas using more microscopic approaches 
such as density-functional theory. The density (or magnetisation) profile for such 
systems is of course highly inhomogeneous and even analysis of a Landau-like model 
would require numerical methods. Beyond mean-field Ising model simulations could test 
the role of the wetting parameter in determining the scale of the droplet height. Finally 
we hope that at least some of our predictions are open to experimental verification from 
observation of droplet shapes in the micro-metre range. 

E.D.Macdonald acknowledges the EPRSC for a studentship and CR. acknowledges 
the financial assistance of the E.G. under contract ERBFMBICT983229. 
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Appendix A. 

We show here that, in two dimensions, the equation 

^^i = Ce-~^\ (A.l) 

(where C and k are real numbers) is conformally invariant. In particular, we demonstrate 
that Eq. (|A.1|) is invariant under the following conformal transformation 

{x,y) — > iu,v) 

^-k£(x,y) ^ ^-kl{u,v) jy^'ji^b (A. 2) 

where {u,v) are given by the relation (in the complex plane) u + iv = w{z), z = 
X + iy, w{z) is any analytic function, \w'\ is the local length rescaling factor of the 
transformation and Xb = 2. 
If we define 

V = 



J 




substitution of ([A2D into ( |AT]) yields 



_ ^ V' log J = C e-'^^(«''') J^^/\ (A.3) 
2k 

Using analytic properties of the conformal transformation 

du dv du dv 

dx dy' dy dx 

and standard methods, it can be shown that log J = and = J V . Therefore, 
eq. (|A.3|) transforms into 



which, in turn, recovers (|A.1|) for Xf, = 2. 
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Figure 1. Schematic illustration of a heterogeneous planar substrate containing a 
finite-size domain of substrate type A which is completely wet by liquid at bulk two- 
phase coexistence corresponding to local contact angle 9a = 0. Outside this domain 
the substrate is of type B and is only partially wet so 6b > 0. All lengthscales are 
considered larger than the bulk correlation length. 
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Figure 2. Heterogenous substrate made from the union of two semi-infinite surfaces 
of type A and B. The mean interface height grows to infinity as we move further into 
the completely wet domain. For short ranged forces the function e~'^^^'> has the scahng 
form e"**^^^ o<:y~^ far from the y = line. 
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Figure 3. Schematic illustration of a droplet cross-section on a heterogenous striped 
domain of finite- width L. For short-ranged forces the analytic expression for the scaling 
of the droplet shape is obtained by conformally mapping the semi-infinite substrate 
shown in Fig.^ 
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Figure 4. Comparison of tiie scaled interfacial lieiglit profiles for the strip and 
disc geometries for systems with dispesion forces. The disc radius R is equal to the 
strip half-width. The height is measured in appropriate dimensionless units whilst the 
scaling variable x/L denotes both the scaled distance across the strip and the radial 
distance across the disc. 
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Figure 5. Section of an interfacial configuration occurring in a wedge geometry. As 
the critical filling transition temperature Tw is approached the mean mid-point height 

roughness and anisotropic correlation lengths ^2;, diverge continuously. The 
extreme anisotropy of fluctuations means that the transition is quasi-one dimensional. 
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Figure 6. Schematic illustration of an intcrfacial configuration in the wedge geometry 
viewed on a scale larger than Fluctuations in the local height across the wedge 
have disappeared leaving only the quasi-one dimensional undulations of the mid-point 
height along the system. 




Figure 7. A heterogenous wedge containing a finite-size region A which satisfies the 
local filling condition dA{TF)=a. Outside this strip, in the region B the wedge is only 
partially filled. 



